Anisotropic spatially heterogeneous dynamics in a model glass-forming binary mixture 
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We calculated a four-point correlation function G4,(k,f;t) and the corresponding structure factor 
Si (k, q; t) for a model glass-forming binary mixture. These functions measure the spatial correlations 
of the relaxation of different particles. We found that these four-point functions are anisotropic and 
depend on the angle between vectors k and f (or q). The anisotropy is the strongest for times 
somewhat longer than the (3 relaxation time but it is quite pronounced even for times comparable 
to the a relaxation time, r a . At the lowest temperatures 5*4 (k, q; r Q ) is strongly anisotropic even for 
the smallest wavevector q accessible in our simulation. 

PACS numbers: 64.70.Pf, 61.20.Lc 



The hypothesis that there is a growing dynamical cor- 
relation length which accompanies the glass transition, 
and that this correlation length is associated with the 
dramatic slowing down of a supercooled liquids' dynam- 
ics, has recently prompted many computational and the- 
oretical PjSSSIMIaiHSUj and some experimental 
[lfl IrH ] investigations. Several studies have been moti- 
vated by the observation that the relaxation of the su- 
percooled liquid involves the correlated motion of clus- 
ters of particles, and the size of the clusters increases 
with decreasing temperature Since two-point cor- 

relation functions, for example the van-Hove correlation 
function, do not provide any information about the cor- 
related motion of particles, four-point correlation func- 
tions were introduced. The dynamic correlation length 
was determined by studying the spatial decay of these 
four-point correlation functions (or the small wavevec- 
tor dependence of their Fourier transforms), and several 
studies have found that this dynamic correlation length 
increases with decreasing temperature 0, 0, IE EH • 

The four-point correlation functions are usually as- 
sumed to be isotropic (in some cases they are isotropic 
by construction Q). However, several researchers have 
noted that on the (3 relaxation time scale the correlated 
motion of particles is not isotropic. Doliwa and Hcuer 
|l4j observed pronounced anisotropy in the f3 relaxation 
regime for a hard sphere system. Glotzer et al. reported 
that the motion of " mobile" particles for a model binary 
Lennard-Jones liquid is not isotropic, but the particles 
move in patterns which they referred to as "string-like" 
0,01- In other words, a particle would follow another 
particle in the liquid in quasi-one-dimcnsional "strings." 
These observations suggest that the correlated motion of 
the liquid's molecules is anisotropic during (3 relaxation. 
Therefore, an appropriately defined four-point correla- 
tion function used to measure this motion should also be 
anisotropic, at least in the f3 relaxation regime. In this 
Letter we discuss a four-point correlation function and 
the corresponding structure factor which are anisotropic 
on the time scale of both j3 and a relaxation. 

To calculate these correlation functions we performed 



Brownian dynamics simulations of the Lcnnard- Jones 
80:20 binary mixture, which was first introduced by Kob 
and Andersen Il7i. The details of the simulations are de- 
scribed in Refs. |l8lll9j |. In this Letter we only present re- 
sults for the larger and more abundant A particles. Thus, 
all sums over particles in the formulae below run over the 
A particles only. In the figures we present the distance 
as r/(TAA, wavevector dependence as qaAA, and time as 
tDg/a AA ; & aa is the Lennard-Jones radius for the inter- 
action among the A particles and Dq is the short time 
diffusion coefficient, which is temperature dependent in 
our simulations. 

A four-point correlation function that we study char- 
acterizes correlations between relaxation of different par- 
ticles. Consider the function 



# n (fc;t)= e *P«( t )-^(°)] ) 



(1) 



where f n (t) is the position of particle n at a time t. The 
ensemble average of F n (k; t) is equal to the self intermedi- 
ate scattering function F s (k;t). The four-point function 
G4 measures correlations between the microscopic self- 
intermediate functions pertaining to different particles 
separated at t = by vector r (see Fig. QJ, 



F n (k;t)F m {-k; t)S[f - f nm (0)} 



G 4 (k,r;t) = ^Y, 

(2) 

where r nm = f n — r m , V is the volume, and TV is the 
number of particles. In this work we fixed \k\ — 7.25, 
which is around the value of the first peak in the static 
structure factor [ljj, |2fJ. Note that G±(k, r; 0) = g(r) 
where g(r) is the pair correlation function. 

The microscopic self-intermediate functions F n (k;t) 
are sensitive to particles' motions along the wavevector 
k. Thus, the four-point function Gi(k, r; t) measures cor- 
relations between particles' motions along k. We expect 
that for small values of \r | the anisotropy of these corre- 
lations is most pronounced for r parallel to k. 

To investigate local correlations on the a relaxation 
time scale |2jJ we calculated projections of Gi(k, r;t) 
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FIG. 1: G 4 (k,f;t) is a pair correlation function in which 
contributions of the individual particles are weighted by their 
microscopic self-intermediate functions F n (k;t). 



onto the Legendre polynomials, 

L n (k,r;t) = f G 4 (k,f,t)P n (tf)d?, (3) 



where P n is the n th Legendre polynomial, k = k/k, f = 
f/r, and dr denotes integration over a unit sphere. If 
G 4 (k,r;t) does not depend on the angle between k and 
f, then L n (k,r;t) will be zero for n > 0. Note that 
the imaginary part of G 4 (k,r;t) is not zero, thus there 
are non-zero real and imaginary parts to L n (k,r;t); by 
symmetry, the imaginary part of L n (k, r; t) is identically 
zero for even n, and the real part is zero for odd n. 

In Fig. |21 we present results for L (k,r;r a ), Im 
.Li (k, r\ Tr,), and L2 (k, r; r Q ) where r a is the a relaxation 
time |21|. As expected, even on this long time scale there 
are strong local correlations between particles' relax- 
ation. The height of the first peak of the isotropic compo- 
nent of G4, Lo(k,r;r a ), increases with decreasing tem- 
perature. Furthermore, we find pronounced anisotropy 
of the local correlations. The correlations revealed by 
G 4 (k,f;T a ) are the strongest when vectors k and r 
are parallel. The amplitude of the anisotropic part of 
G 4 (k,r;T a ) is, roughly speaking, temperature indepen- 
dent for T < 1.0. 

To examine correlations between particles' relaxation 
for large distances we turn to the structure factor corre- 
sponding to the four-point correlation function. We de- 
fine Si{k,q\t) = l + (N/V)H 4 (k,q;t), where H 4 (k, q\t) is 
the Fourier transform of G 4 (k, r; t) — F^(k; t). For q ^ 
we have 



S 4 (k,q;t) = — ^ 



F n (k;t)F m (-k;tyf- f 



(4) 



Four-point wavevector-dependent functions similar to 
S 4 {k,q;t), usually evaluated at t = r Q , have been 
used t 

QUI 

First, in Fig. [3] we examine projections of S 4 (k,q;t) 



ireviously to determine dynamic correlation lengths 
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FIG. 2: The projections Lo(k,r;T a ), Im Li(k,r;r a ), and 
L 2 (k,r;T a ) for T = 0.45, 0.55, and 1.00. 



onto the Legendre polynomials, 

In + 1 f - ^ - * 

I n (k,q;t)= 4 / S 4 (k,q,t)P n (k ■ q)dq, (5) 

at the a relaxation time. We do not show Ii(k, q; t) since 
it is small for qoAA < 5. Next, in Fig. 01 we examine 
Si(k, q; r a ) for different angles 9 between k and <£In ac- 
cordance with numerous earlier investigations [ij, |2j, |3, |^] , 
we find that the large long-range correlations between 
particles' relaxation increase with decreasing tempera- 
ture. Surprisingly, we find that even the long-range 
correlations on the a relaxation time scale are strongly 
anisotropic. At low q the correlations are strongest for 
q _L k. In direct space this corresponds to strong cor- 
relations for particles separated by a vector r which is 
parallel to k. 

The anisotropy of S 4 (k,q;t) makes the determination 
of a single dynamic correlation length difficult. We tried 
to fit the low q values of S 4 (k, q; t) for a fixed angle 9 be- 
tween k and q to several different functional forms, but 
could not find a form where the results were reasonable 
for T < 0.50. Specifically, fits to an Ornstein-Zernickc 
formula were very poor. The correlation length deter- 
mined from the fits was either larger than half the simu- 
lation cell, or we had to fix the unknown q = value of 
S 4 (k, q; t). Larger systems need to be simulated in order 
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FIG. 3: The projections Io(k,q;r a ), and l2(k,q\T a ) for T = 
0.45, 0.55, and 1.00. 



to better understand the low q behavior of Si(k, <f; t). 

While it has been previously shown that heterogeneous 
dynamics on the (3 relaxation time scale is anisotropic 
. the strong anisotropy of Si(k, q\T a ) was a surprise. 
To investigate the time dependence of the anisotropy 
we examined the ratio S±{k, q±; t)/Si{k, qii ; t), where 
S4(k,q±;t) is calculated for the smallest wavevector al- 
lowed due to the periodic boundary conditions (\q±\ = 
2-k/L) perpendicular to k and Si(k, q« ; t) is calculated for 
the smallest wavevector (\q^ \ = 2-k/L) parallel to k. This 
ratio is shown in Fig. for T = 1.0, 0.55, and 0.45. The 
arrows in the figure indicate the a relaxation time and 
the f3 relaxation time, rg. The [3 relaxation time was de- 
termined by finding the first inflection point of \n[F s (k; t)\ 
as a function of t. For T = 1.0, this inflection point does 
not exist. The peak in the ratio S^k, q±;t) / S^k, q\\;t) 
(i.e. the maximum anisotropy) occurs between the (3 and 
the a relaxation time, but is closer to the /3 relaxation 
time for the lower temperatures. However, the correlated 
motion is still strongly anisotropic around T a . 

Finally, we compared Si(k,qo',t) for the smallest 
wavevector allowed due to the periodic boundary con- 
ditions, <fo| = 27r/L, with the four-point susceptibility 



F n (k;t)) (F m (-k;t) 



(6) 



As observed before 20], while Si(k,q = 0;t) = Xi(k;t), 
the susceptibility Xi(k',t) is ensemble-dependent, and in 
a constant N ensemble lim^o Si(k, q)t) ^ Xi(k',t). In 
Fig. we show the correlation function Si(k, qo]t) and 
the susceptibility ^(fcjt) for several angles 8 between 
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FIG. 4: The correlation function Si{k,q\T a ) for T = 0.45, 
0.55, and 1.0. The angle 9 is the angle between k and q. 
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FIG. 5: The ratio S4k,q±;t)/Si(k,qf,t) for T = 1.0, 0.55, 
and 0.45. The arrows indicate the a and the /3 relaxation 
times. 



vectors <fo and k for T = 1.0, 0.55, and 0.45. For a fixed 
9 there is a peak in S4,(k, qo',t) around r a . The exact 
position of the peak depends on the angle between <fo 
and fc, and occurs at a time larger than r a for qo parallel 
to k, but smaller than the r a for <fo perpendicular to k. 
It is clear from Fig. that the time dependence of the 
correlation function Si(k, qo] t) and that of the suscepti- 
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FIG. 6: The four-point correlation function S4,(k, qo;t) as a 
function of time for T=0.45, 0.55, and 1.0. The solid lines are 
for the angles 0°, 30°, 45°, 60°, and 90° listed from bottom 
to top. The dashed lines represent the angles 120°, 135°, 
150°, and 180° listed from top to bottom. The dark dashed 
dotted line is Xi(k',t) and the vertical dotted line marks the 
a relaxation time. 



bility xt(k',t) are quite different. This is significant be- 
cause time dependence predicted by recent extensions of 
the mode-coupling theory |23j for the correlation function 
S±(k,q\t) are usually checked a gain st simulation results 
for the susceptibility X4(fc;i) |&|24|. 

In summary, we have shown that there is a pronounced 
anisotropy of the correlations of particles' relaxation in 
a model supercooled liquid even on the time scale of the 
a relaxation time. This anisotropy will need to be ad- 
dressed in theoretical descriptions of heterogeneous dy- 



namics of supercooled liquids. In particular theoretical 
input is needed in order to elucidate the relation be- 
tween anisotropic four-point correlation functions and a 
dynamic correlation length. Finally, we expect that the 
anisotropy of the four-point correlation functions exam- 
ined in this Letter will help to differentiate between com- 
peting theoretical descriptions of heterogeneous dynam- 
ics. 

We gratefully acknowledge the support of NSF Grant 
No. CHE 0517709. 
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